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A network of N elements is studied in terms of a deterministic globally coupled map which 
can be chaotic. There exists a range of values for the parameters of th e map where the number of 
different macroscopic configurations jV(JV) is very large, Af(N) ~ exp yj c(a)N and there is violation 
of self averaging. The time averages of functions, which depend on a single element, computed over a 
time T, have probability distributions that do not collapse to a delta function, for increasing T and 
JV. This happens for both chaotic and regular motion, i.e. positive or negative Lyapunov exponent. 
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Glassy behavior is still one of the central subjects of 
modern physics Q. Among the most interesting fea- 
tures are the violation of time traslation invariance, and 
fluctuation-dissipation theorem, even for very long times 
||. This follows from the breaking of ergodicity in the 
phase space caused by a rough energy landscape. In spin 
glasses, the paradigmatic system, the complex dynamics 
follows from frustration and randomness. It is, however, 
obvious that these are not a prerequisite for glassy be- 
havior, think for example to ordinary window glass. Re- 
cently an increasing amount of work has been devoted 
to the study of glassy behavior in non random systems 
||. In this case the complex behavior follows from the 
existence of a large number of metastable states which 
slows down the dynamics. 

We shall concentrate on systems of globally coupled 
elements without any intrinsic disorder. Dynamical sys- 
tems of this type have been studied recently, in great 
details, in many fields 0, and used to model a wide class 
of phenomena, such as Josephson-junction array, evolu- 
tion dynamic of ecosystems, neural-networks for neuro- 
dynamics. 

The interesting point is that the macroscopic descrip- 
tion of these systems may be extremely rich. The variety 
of the macroscopic description can be taken as an indi- 
cation of complexity || . 

In this Letter we shall study a network of elements - 
a globally coupled map - given by: 

JV 

Xi (n + 1) = (1 - e)f( Xi (n)) + ^ ]T /fo(n)) (1) 

3=1 

where n is the discrete time variable and i = 1, • • • , JV is 
the index identifying the elements. For f(x) we choose 
the logistic map 



a violation of the law of large numbers: a subtle coher- 
ence among the elements Xi induces anomalous fluctua- 
tions of the mean field 



1 — ax 



(2) 



The system ([T]), that has been thoroughly studied by 
Kaneko p], displays a very rich behavior. Interestingly, 
for certain values of the parameters a and e, one observes 



1 N 



(3) 



so that the mean square deviation, (h 2 ) — (h) 2 , decays as 
JV _/3 , with (3 < 1. Moreover, the elements split into clus- 
ters whose probability distributions resembles, for suit- 
able values of the parameter, to those derived from the 
theory of spin glasses [Q . 

In this letter we show that the analogy of this sys- 
tem with spin glasses, previously proposed by Kaneko, 
is very deep: for a certain range of values of the control 
parameter a, the system exhibits a very large numbers of 
"macroscopic" states and a violation of selfaveraging, in- 
dicating that the system is "complex" in the above stated 
sense. 

Because of the global coupling in ([!]), one can have the 
phenomenon of clusterization: if, at a certain time n, the 
variables of some elements take on identical values, they 
will maintain equal values for every n > n. Therefore, 
it is possible to introduce a "macroscopic" description 
of the system in terms of the number JVj of elements in 
the i th cluster and to identify a macroscopic state by the 
set of integers {JVj} = [Ni, JV 2 , . . . , JV m | J2ti W = 
that characterize the population of each cluster. All the 
microscopic states with the same cluster structure, i.e. 
with identical values for the {Ni}, but differing either by 
a permutation of some Xj or by the value of the variables 
Xj in a particular cluster, belong to the same macroscopic 
state (by definition). 

We start by estimating the number of relevant macro- 
scopic states of the system as a function of JV, for differ- 
ent values of a. This is done by generating for each value 
of JV and a a set M. 1 of initial conditions x^ a '(0), 
with a = 1,2, ■■■,M, selecting at random each com- 
ponent x^ (0) with uniform distribution in the interval 
[—!,+!]; we used M = 10 4 -f5x 10 4 . The numerical 
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results reported here were obtained with e = 0.1. For 
each initial condition the map (|l|) is iterated and the fi- 
nal "cluster structure" recorded. In all cases we have 
discarded 6 • 10 4 thermalization steps before performing 
measurements. From each set of initial conditions we ob- 
tained a "cluster structure" , i.e. a set j-^j"^ j > and, from 
the frequency of occurrence of every set, we estimated the 
probability, P {N\, N2, ■ ■ . , N m ) , to observe the macro- 
scopic state (Ni, N2, ■ . . , N m ). To take into account only 
the relevant states, i.e. to neglect the configurations with 
very small probability, we compute the entropy 

S(N) = - Pn(N 1 ,N 2 ,... 1 N„ 1 ) x 

(iVl,JV 2 ,...,iV m ) 

lnP^A^A^...,^) (4) 

from which we extract the number of macroscopic con- 
figurations with non negligible probability 

Af(N) = e s(N \ (5) 

A clusterization of N elements in the language of number 
theory is a partition of the integer N, i.e., a possible way 
to obtain TV by summing integers - not necessarily all 
different. An upper bound for S(N) and a guess for how 
it scales with the system size can be obtained noting that 
the asymptotic estimate for the number of partitions of 
iV,A/- p (A0,is § 

tt^2A73 

From this we may guess that 

S(N) 2 ~ c(a)N + cost (7) 

with c(a) < c max = (2/3)7r 2 ~ 6.58. The constant c(a) 
may be thought as an indicator of the fraction of the "ac- 
tive" or "relevant" elements in the system. In fact, if one 
defines the number N act by the relation 

Af(N) ~ e*>/ 2N <"*/ 3 , (8) 
it readily follows that 

c(a) = Nact_ , g s 

Cmax N 

In Fig. @ we display S(N) 2 as a function of N for dif- 
ferent values of a. For each initial condition we checked 
that the cluster structure, found at the end of the ther- 
malization, is a "permanent" structure, that was exactly 
the same after 6 x 10 4 and 1.2 x 10 5 successive iterations. 

Figure shows c and the mean value (A) of the maxi- 
mum Lyapunov exponent of the system averaged over the 
initial conditions. The bars on (A) indicate the variance 
<7at. For a > 1.63 the variance ctn goes to zero as N in- 
creases and the probability distribution of the Lyapunov 



exponents, Pjy(A), approaches a delta function, around 
some A > 0, as N — > 00. On the other hand, for a < 1.63 
and for N — > 00 ct^v approaches a non zero value, while 
P/v(A) approaches a nontrivial asymptotic distribution. 

It is evident that for a in the range 1.5 -7- 1.7, the 
system has a very large number of macroscopic states, 
while outside this range one detects only very few macro- 
scopic states. We stress that we observe many macro- 
scopic states - both chaotic and regular - not only many 
microscopic states, as in Ref. || where, in a system of 
N globally coupled nonlinear oscillators, (N — 1)! stable 
limit cycles are found. Moreover, the maximum of c(a) 
is attained in the region where (A) changes sign. 

The above results and the observation of the his- 
tograms of the maximum Lyapunov exponent, suggest 
the existence of four distinct regions in the behavior 
of the system. These are roughly equivalent to those 
identified by Kaneko. We stress, however, that whereas 
Kaneko's classification is based on microscopic considera- 
tions, ours follows from the observation of a macroscopic 
behavior. The four regions are: 

1. a < 1.5, standard regular motion, with (A) < and 
few macroscopic states. 

2. 1.5 < a < 1.63, "glassy" regular motion, i.e.: 
(A) < and many macroscopic states. 

3. 1.63 < a < 1.70, "glassy" chaotic motion, with 
(A) > and many macroscopic states. 

4. a > 1.70, standard chaotic motion, with (A) > 
and few macroscopic states. 

The analogy with the glassy behavior in disordered sys- 
tems is not limited to the number of macroscopic states: 
in fact, in its glassy phase, system (Q) is characterized 
by a violation of the selfaveraging. Given an initial con- 
dition, x(0), we may compute the time average, over a 
time T, of the variable xf. 

1 T 

n=l 

or of any other function of a single element. In a stan- 
dard situation, the probability distribution of {x%)t is 
expected to collapse to a delta function as T and N in- 
crease, that is, we expect that the time average over a 
very long time record of one element chosen at random 
inside a large system gives a complete description of the 
one- variable properties of the system. 

Figure 3a shows the probability distribution of (xi)T, 
obtained by gathering the average values for different i 
and different initial conditions, in the phase with stan- 
dard chaos. By increasing the observation time T, the 
variance decreases and one recover the usual behavior. 
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Figure 3b shows the situation in the case of the chaotic 
glassy phase: the increasing of T does not modify the 
nontrivial structure of the distribution. A similar behav- 
ior is observed in the regular glassy phase. 

These results, which are size independent, show that, 
in the glassy phase, the ergodicity is broken and the self- 
averaging property does not hold, in other words: the 
time average on a single element is not sufficient for a 
complete characterization of the one-element statistics. 

We conclude with a brief discussion about the role of 
noise. Of course, for finite values of N, the presence 
of noise will restore ergodicity after a very long time, 
which depends on N. Therefore it is interesting to con- 
sider the behavior of the system subject to noise on large 
but finite times. In Ref. || for large N even a small 
amount of noise is sufficient to induce transitions among 
the different attractors. On the contrary we present a 
different picture. Within our approach, the scenario de- 
scribed above for the system ([!]) is very robust with re- 
spect to noise. In fact, the inclusion of a noise term 
arji (n) - where T)i are independent random variables uni- 
formly distributed in the interval [—1, 1] - to the r.h.s. 
of equation (Q) does not change the probability distribu- 
tions of {xi)T showed in Fig. 3a for noises strength as 
large as a » 10~ 2 . Therefore the macroscopic states are 
stable for long times. This feature closely resemble to 
the behavior observed in glassy systems below the glassy 
transition temperature. 

In conclusion we have shown that even simple mod- 
els which do not contain any intrinsic disorder can have 
an extremely rich behavior at a macroscopic level. A 
behavior typical, e.g., of glasses. This scenario is very 
robust with respect to small external noise, indicating a 
"genuine" property of the system. 
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FIG. 1. Squared entropy S 2 as a function of the number of 
globally interacting elements N, for values of the parameter 
a equal to 1.50, 1.55 and 1.63. The lines are the linear best 
fits. 



FIG. 2. Entropy coefficient c (boxes) and average Lya- 
punov exponent (A) (triangles) as functions of a. The values 
of (A) and of crjv are computed with N ranging from 100 up 
to 3200, and with M = 10 4 initial conditions. 
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FIG. 3. Probability distribution of (xi)t, at various obser- 
vation times T in the region of standard chaos - a — 1.80 - (a) 
and in the region of glassy chaos - a = 1.67 - (b). The differ- 
ent curves correspond to averages performed, after a transient 
of 6 x 10 4 steps, on time intervals T = 0.75 x 10 4 , 1.5 x 10 4 , 
3 x 10 4 , 6 x 10 4 , 12 x 10 4 . For each curve M = 10 3 initial 
conditions were considered with N = 300. The histograms do 
not change qualitatively by changing N up to 1200. 
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